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Abstract 

We explore the origins of slow dynamics, aging and glassy rheology in soft and 
living matter. Non-diffusive slow dynamics and aging in materials characterised by 
crowding of the constituents can be explained in terms of structural rearrangement 
or remodelling events that occur within the jammed state. In this context, we in- 
troduce the jamming phase diagram proposed by Liu and Nagel to understand the 
ergodic-nonergodic transition in these systems, and discuss recent theoretical at- 
tempts to explain the unusual, faster-than-exponential dynamical structure factors 
observed in jammed soft materials. We next focus on the anomalous rheology (flow 
and deformation behaviour) ubiquitous in soft matter characterised by metastabil- 
ity and structural disorder, and refer to the Soft Glassy Rheology (SGR) model 
that quantifies the mechanical response of these systems and predicts aging under 
suitable conditions. As part of a survey of experimental work related to these is- 
sues, we present x-ray photon correlation spectroscopy (XPCS) results of the aging 
of laponite clay suspensions following rejuvenation. We conclude by exploring the 
scientific literature for recent theoretical advances in the understanding of these 
models and for experimental investigations aimed at testing their predictions. 
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1 Introduction 

Certain out-of-equilibrium soft materials, including foams [1], concentrated 
emulsions [2] and fractal colloidal gels [3], exhibit slow relaxation processes 
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that often resemble the glassy dynamics observed in disordered hard materi- 
als such as polymer melts and spin glasses. The properties of many of these 
disordered soft systems evolve slowly with time, a behaviour commonly termed 
'aging'. Aqueous foams, for example, age due to coarsening and bubble rear- 
rangement processes [4]. Polymer glasses [5], microgel solutions [6] and aqueous 
clay suspensions [7,8] are some more examples of aging systems. In spite of 
their apparent diversity, these materials are all structurally disordered and 
out-of-equilibrium. The existence of certain universal features in the evolution 
of the mechanical responses of these systems indicates an underlying generic 
mechanism that is responsible for the observed behaviour. To this end, a num- 
ber of theoretical and experimental studies of the evolution of aging systems 
- including calculations of magnetic suceptibilities in frustrated spin systems 
[9], dielectric measurements of memory in molecular glasses [10], and measure- 
ments of the dynamic correlation functions of colloidal gels [3] - have noted 
commonality in the dynamics under study with seemingly disparate systems. 

In addition to being sources of challenging fundamental problems in their 
own right, soft materials can serve as excellent model systems in the under- 
standing of hard condensed matter as they often display several dynamic and 
thermodynamic properties typically associated with the latter. For example, 
recent work on melting of colloidal crystals has elucidated the importance of 
pre-melting at grain boundaries [11]. Another example is the mode coupling 
theory (MCT). Originally applied to hard sphere liquids, MCT accounts for 
many features observed in the colloidal glass transition including the pres- 
ence of a two-step relaxation in the dynamical correlation function [12,13,14]. 
Based on this success, researchers have sought to extend application of the 
theory to the glass transition in molecular glasses and polymers [15,16]. Ow- 
ing to their larger sizes and slower relaxation times when compared to, say, 
atomic fluids, the structure, dynamics and mechanical response of soft mate- 
rials can be studied experimentally by common laboratory techniques such as 
microscopy, light scattering and rheology. The interaction between the indi- 
vidual constituents of soft materials can be tuned easily [17] and this feature 
is exploited for example in the verification of the predictions of the MCT for 
repulsive and attractive colloidal glasses [18,19]. For a more complete discus- 
sion, we would like to refer the readers to a recent review article by Cipelletti 
and Ramos [20]. 

The presence of a glassy response, which is characterized by very slow, history- 
dependent dynamics, implies a loss of ergodicity. The transition of colloidal 
suspensions and granular media from an ergodic to a nonergodic state can be 
understood in terms of the jamming of overcrowded particles, which results in 
kinetic arrest and ultra-slow relaxation processes [21]. For dense materials such 
as glasses, sandpiles and foams, Liu and Nagel [22] suggest a 'jamming' phase 
diagram that considers the possibility of unjamming the system [i.e driving it 
towards ergodicity) by increasing its temperature, by applying appropriately 



2 



large shear stresses, or by decreasing the density of its constituents. Univer- 
sal non-diffusive slow dynamics and aging have been observed in a variety of 
jammed soft materials, for example, in colloidal suspensions [2,3,7,8,23,24], 
concentrated emulsions [2], micellar polycrystals [2] and lamellar gels [2,25]. 
Glassy dynamics also occur in biologically relevant systems. Examples include 
the kinetic arrest and scale-free rheology observed in the cytoskeleton of a liv- 
ing cell [26,27] and the ultraslow relaxation processes in models of the adaptive 
immune response system [28]. 

In the specific case of disordered soft materials that can flow in response to 
macroscopic strains, metastability can result in anomalous rheology that is 
understood in terms of the soft glassy rheology (SGR) theory proposed by 
Sollich et al. [29,30]. For those systems whose thermal energy is not sufficient 
to achieve complete structural relaxation, SGR predicts identical power law 
frequency dependences of the low-frequency elastic and viscous moduli and 
identifies the exponent of this power law with a mean field noise temperature. 
The model predicts a glass transition, with the glass phase being characterised 
by the aging of the viscoelastic moduli and the presence of yield stresses. Many 
of these features have been observed experimentally in concentrated microgel 
solutions [31], foams [32], paint [33], compressed emulsions [34], a thermotropic 
liquid crystal with quenched disorder [35] and in the cytoskeleton of living cells 
[27]. 

In the remainder of this article, we discuss some theoretical models that de- 
scribe the slow dynamics, aging and glassy rheology in soft materials and cite 
experimental results that verify the predictions of these models. In doing so, 
we have organised the paper in the following manner: section 2 contains an 
overview of the jamming phase diagram proposed by Liu and Nagel [22] fol- 
lowed by a brief description of theoretical attempts to explain the non-diffusive 
slow dynamics and aging in soft materials in the jammed state. We next focus 
on the SGR theory [29,30] that quantifies the flow and deformation behaviour 
of soft glassy materials (SGMs). Section 3 reviews experimental work on the 
slow dynamics and aging in soft and living matter, followed by observations 
of anomalous rheology in many of these systems, with some emphasis on our 
recent work on the aging behaviour in synthetic colloidal laponite clay sus- 
pensions [24] and the observation of soft glassy rheology in a thermotropic 
smectic liquid crystal confined within the pores of a compliant colloidal gel 
[35] . This is followed by a brief summary and concluding remarks in section 4. 



2 Models 



In this section, we present brief overviews of some models that are particularly 
relevant in the description of the slow dynamics, aging and glassy rheology in 
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soft and living matter. These include the jamming picture, the prediction of 
ultraslow non-diffusive relaxation processes in aging soft matter and the SGR 
model. Experiments that verify the predictions of these theories are summed 
up in section 3. 



2.1 The jamming phase diagram and the compressed exponential dynamical 
structure factor in some jammed systems 

An approach to understanding the ergodic-nonergodic transition is provided 
by the concept of jamming [21]. Jamming of overcrowded particles results in 
kinetic arrest and signals a fluid-solid transition. The jamming phase diagram 
provides a unified picture for the ergodic-nonergodic transition in systems as 
diverse as molecular glasses, foams, granular matter and colloidal suspensions 
[22] . If a system arrests as soon as it is able to support an external load, such 
that it responds elastically to some compatible loads while undergoing plastic 
rearrangments under incompatible loads, it is classified as 'fragile' [36,37,38]. 
Jamming of fragile matter leads to slow dynamics because of the kinetic con- 
straints it imposes on the motion of the individual constituents. Cates et al. 
draw a distinction between fragile materials and conventional solids by noting 
that in the case of the former, applied stresses are supported by force chains 
that can collapse under the action of an incompatible load [38,39]. This can 
lead to unjamming. The jamming picture has been subsequently expanded by 
Liu and Nagel [22] to include jammed systems such as foams and granular 
matter which are characterised by yield stresses and are therefore not fragile 
according to the definition of Cates et al. [38]. In addition to fluidising a sys- 
tem by the application of shear stresses as suggested in [38], the authors in [22] 
suggest the possibility of unjamming a jammed system by increasing its tem- 
perature (in the case of, say, molecular glasses) or by decreasing the density of 
its constituents (in the case of, say, aqueous foams). In this way, they suggest 
a natural connection between the glass transition in supercooled liquids and 
the fluid-solid transition in athermal systems such as foams and granular mat- 
ter. To encompass these ingredients, the authors introduce a phase diagram 
(Fig. la) for jammed systems whose axes depend on temperature, load and 
the inverse of the density, such that appropriate changes in these parameters 
can result in the onset of the unjamming transition. 

Ultra-slow dynamics and aging observed in colloidal gels [3] , clay suspensions 
[7,8,24], compressed emulsions [2] and multilamellar vesicles [25] are under- 
stood in terms of the relaxation of internal stresses that are built into the 
sample at the jamming transition. The dynamical structure factors (DSFs) 
characterising these systems are strongly dependent on sample age (quanti- 
fied in terms of the waiting time t a since sample preparation or the onset of 
jamming) and exhibit non-diffusive, faster-than-exponential forms. The relax- 
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Fig. 1. The jamming phase diagram proposed in [22] is depicted in a). The jamming 
phase diagram for attractive colloidal particles constructed on the basis of exper- 
iments reported in [59] is depicted in b), where <j) is the colloidal volume fraction, 
U is the interparticle attractive energy whose scale is set by the thermal energy 
/cgT and a/a Q is a suitably normalised applied stress. The shaded regions in both 
diagrams indicate jamming. The phase boundaries are qualitative, and may vary 
between systems. Figures adapted from [22] and [59]. 

ation time t s characterising these slow dynamics scales with wave vector q 
as r s ~ q^ 1 . In the case of a fractal colloidal gel which undergoes syneresis 
(microcollapses in its structure) [2,3], the relaxation process that gives rise to 
the observed decay of the DSF is modeled in terms of the elastic response to 
local stress sources that appear at random in the medium. Each stress source 
is approximated as a force dipole, such that for a small number density of 
stress sources distributed randomly in space, calculations predict the relation 
r s ~ q^ 1 and a compressed exponential form for the DSF: f(q, t) ~ exp[— (qt) p ] 
with a compressing exponent p ~ 1.5. These trends imply ballistic strain in 
the sample characterised by a velocity distribution function W(v) 
response to the stress [2]. Allowing for a finite collapse time 9 for each micro- 
collapse event and writing a central force elasticity equation for the strain field 
around a force dipole, Bouchaud and Pitard [40] calculated several regimes for 
the DSF, with the compressed exponential behaviour f(q,t) ~ exp[— (gt) 15 ] 
occuring only at earlier times. With a coupling between stress dipoles incor- 
porated into the model, the growth of t s with t a is also computed, giving a 
t s ~ t^/ 3 dependence in this early time regime. The compressed exponential 
form of the DSF and the inverse linear dependence of the relaxation time on 
the wave vector are seen in experiments in a wide variety of soft materials. A 
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sublinear dependence of r s on t a is also often observed, though the exponent 
of the power law can deviate significantly from the theoretical prediction in 
[40] . This variability indicates a strong dependence of the aging process on the 
structural and dynamical details of the sample under investigation. 

2.2 Trap models and Soft Glassy rheology (SGR) 

Two-step relaxation processes, such as those observed in supercooled liquids 
[15,16], dense colloidal suspensions [13] and concentrated emulsions [41], can 
be understood in terms of the trap model. This model describes the glassy 
dynamics of a particle (or an 'element') trapped in a cage, i.e. in a potential 
well created by the constraints imposed by neighbouring particles, such that 
escape from a cage is possible only by some activated process (for example, by 
thermal activation in the case of supercooled liquids). These potential energy 
minima represent metastable states that are surrounded by high energy barri- 
ers whose heights, in the case of supercooled liquids, increase with decreasing 
temperature. At short times, the particles diffuse in their cages while at longer 
times they escape through activation. The separation of time scales for these 
processes leads to a two-step relaxation. In the glassy state, the escape pro- 
cess becomes inaccessibly slow leading to nonergodic, partial relaxation and 
kinetic arrest. In polyelectrolyte microgel pastes, for example, individual mi- 
crogel beads are trapped in cages formed by the neighboring particles. The 
intensity autocorrelation functions measured for suspensions of these parti- 
cles, which can execute sub-diffusive motion in their cages at short times and 
are trapped (kinetically arrested) at long times, display an incomplete decay 
[42] . Bouchaud used the trap model to understand the relation between aging 
and weak ergodicity breaking in spin glasses [43], where an effective particle 
is trapped in a one-dimensional potential energy landscape, and which, when 
activated, can hop into another trap. Detailed calculations by Monthus and 
Bouchaud [44] show a dynamical phase transition between an ergodic (liq- 
uid) phase at high temperature and a non-ergodic aging (glass) phase at low 
temperature when the density of energy barriers decays exponentially. 

Sollich et al. [29,30] generalised the trap model to understand the low-frequency 
rheology of SGMs. The SGR model proposed by Sollich et al. [29,30] describes 
the flow and deformation behaviour of soft materials such as colloidal pastes, 
foams, emulsions and slurries that are characterised by metastability and dis- 
order. The 'glassiness' of these materials is a consequence of their structural 
disorder and metastability, while their 'softness' arises from their ability to 
flow upon the application of macroscopic strains. SGR [29,30] uses a gener- 
alised trap model to describe SGMs, such that each individual element, which 
moves in a one-dimensional piecewise quadratic potential (Fig. 2), can un- 
dergo 'activated yielding' processes due to structural rearrangement events in 



6 



Ay 

I 



Fig. 2. The potential well picture that forms the basis of the SGR theory for 
metastable and disordered soft systems. The potential energy landscape is charac- 
terised by a distribution of yield energies E and local strains I. As the system is 
sheared beyond a yield point (described by a finite value for the yield strain l y ), 
stress relaxation occurs by the rearrangement of particles (or 'elements') and the 
system settles to a new local equilibrium configuration. The different configurations 
of the system are shown as circles numbered from 1 to 5. A7 corresponds to dis- 
placements when the strain is less than the yield value necessary for the system to 
hop into the next potential energy minimum. When the yield strain is large enough, 
the particle can undergo noise-induced hops into the next minimum, and the en- 
ergy dissipated in such activated yielding processes (denoted by curved arrows) is 
indicated by the solid vertical bars. Figure adapted from [30]. 



the system. Such structural rearrangement events are strongly coupled and 
the resulting interactions are quantified in terms of an effective noise tem- 
perature x. Exact rheological constitutive equations are derived within the 
framework of the SGR model, with the dynamic shear modulus showing a 
power law frequency dependence at low frequencies; for an applied oscilla- 
tory strain j(t) = •ycosut, the elastic modulus G' ~ oj x ~ x for 1 < x < 2, 
while the viscous modulus G" ~ uj x ~ l for 1 < x < 3. The viscous and elastic 
moduli become increasingly frequency independent as x — > 1, with the ratio 
approaching (x — 1) as the glass transition (identified at x = x g = 1) is 
approached. For x < x g , no equilibrium state exists, implying weak ergod- 
icity breaking and the presence of aging phenomena. The flow curve is then 
characterised by an yield stress a y and exhibits the Herschel-Bulkley form: 
a — o y ~ 7 1_a; . This model predicts aging at x < 1 as long as a y is not ex- 
ceeded. For appropriately high values of the noise temperature x, Maxwellian 
relaxation behaviour (characterised by G' ~ u 2 at x > 3 and G" ~ u> for x > 
2) and a power law flow curve (a ~ 7 i '~ 1 for x > 1) are recovered. 
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There have been numerous other attempts to describe the rheology of SGMs 
using such theoretical approaches. For example, Berthier formulated a mi- 
croscopic non-equilibrium MCT that describes the yield stresses, flow hetero- 
geneities and activated processes in an SGM that is characterised by a strong 
coupling between its slow dynamics and the shearing forces [45]. The nonlinear 
rheological behaviour of SGMs characterised by a power law flow curve and 
the existence of a finite yield stress has also been treated within the frame- 
work of the mode coupling model by Hebraud and Lequeux [46]. Lequeux 
and Ajdari [47] calculate a Vogel-Fulcher type divergence for the viscosity in 
their model of SGMs. Calculations by Head et al. indicate that when the noise 
temperature x of the SGM is a decreasing function of the global stress a, the 
flow curve is non-monotonic and the model displays hysteresis [48]. In the 
presence of diffusion, calculations based on hopping models for glassy systems 
[49] show strong violations of the fluctuation-dissipation theorem (FDT) when 
applied to the self-diffusion of small probes in the SGM, while agreement with 
the FDT is restored when the probes are large enough that the surrounding 
medium can be viewed as a continuum. It should be noted that these models 
do not consider the spatial and temporal correlations of the shear-induced 
rearrangement events that are likely to be present in real systems such as 
dense emulsions [50]. The FDT is found to be violated for aging and driven 
systems. However, using a non-equilibrium version of the FD relation, Sollich 
and Fielding show the equivalence between a stationary, weakly driven system 
and an infinitely aging one [51]. 



3 Soft glassy materials: review of experiments 



In this section, we discuss experimental studies related to the predictions of 
the models introduced in the previous section. Most of these experimental 
results are obtained using one of two experimental techniques, viz. rheology 
and photon correlation spectroscopy (PCS), the principles of which we will 
first outline in the following subsection. 



3.1 Experimental techniques 



PCS experiments [52] in the single scattering limit involve the measurement 
of the wave vector dependent intensity autocorrelation function g2(q,t) = 
<I ^ 9 <i(qt')>'^ >T ' ^ere I(q,t') is the scattered intensity measured at time t', 
t is a delay time and T is the total averaging time. The measured g2{q,t) is 
related to the DSF f(q,t) via the Siegert relation g^iq^t) = 1 + P\f(q,t)\ 2 
where f3, the coherence factor, depends on the experimental optics and f(q,t) 
contains information about the dynamics of the scatterers. For example, for a 
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dilute suspension of monodisperse colloidal beads with purely diffusive dynam- 
ics, f(q,t) ~ exp(^) , where r is the characteristic time scale for scatterer 
diffusion and is related to the self-diffusion coefficient D Q according to the 
formula r = {D q 2 )~ l . X-ray photon correlation spectroscopy (XPCS) uses 
synchrotron radiation to probe very high values of q that are inaccessible by 
PCS methods involving visible light [53,54]. For example, XPCS has been 
employed in [24] to study the particle scale dynamics in an aging clay suspen- 
sion, where each particle is a disc approximately f 5 nm in radius and f nm in 
thickness. The dynamics of multiply scattering samples can be accessed with 
diffusing wave spectroscopy (DWS) [55]. Multispeckle DWS measurements of 
probe particles diffusing in aging clay suspensions are used to characterise 
the aging behaviour of the medium and have been reported in [7]. Conven- 
tional PCS experiments use a point detector and therefore require considerable 
time-averaging of the intensity autocorrelation function to ensure reasonable 
statistics. In order to reduce total experimental time, multi-element detectors, 
such as charge-coupled device (CCD) cameras, are used to autocorrelate in- 
tensities at every camera pixel. For isotropic, nonergodic materials in which 
a time average and ensemble average are not equivalent, such multispeckle 
spectroscopy [56] provides a key advantage in that the ensemble average may 
be obtained by an average over pixels corresponding to the same wave vector. 
For these reasons, multispeckle spectroscopy has proved to be an excellent tool 
to probe materials with slow and spatially heterogeneous dynamics [2,3,7,24]. 

In contrast to PCS which measures microscopic properties, rheological exper- 
iments quantify the bulk flow and deformation behaviour of materials. Bulk 
rheology involves the measurements of the viscoelastic moduli and the non- 
linear flow properties of materials [57]. Depending on the magnitude of the 
applied strain, rheological measurements are classified as linear and nonlinear. 
Linear rheology measurements cause no appreciable change in the microstruc- 
ture of the sheared sample, and for a system in equilibrium, the resulting 
linear response functions can be related to dynamic correlation functions for 
the system through fluctuation-dissipation relations. For linear measurements 
involving the application of an oscillatory strain with an angular frequency 
uj and a suitably low amplitude 7, measurements provide the resultant stress 
amplitude o{uS) and the angle 5(u) by which the measured stress is phase- 
shifted with respect to the applied strain. The elastic modulus G' and the 
viscous modulus G" are then computed from this data using the relations 
G'{u) = cos5(cu) and G"(u) = sin<5(cj). For a conventional viscoelas- 
tic material such as a Maxwell fluid, bulk linear rheology provides material 
parameters such as the relaxation time r R and the high frequency plateau 
modulus G [58]. Large values of applied strains, by contrast, can result in 
significant changes in the microscopic structure of the material and can give 
rise to highly nonlinear effects such as shear thinning and thickening, nonzero 
normal stresses etc. [57]. Typical nonlinear rheology experiments include the 
measurement of the flow curve (the stress a vs. the shear rate 7 plot that 
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provides information on the shear rate dependent viscosity 77(7)) and yield 
stresses of materials (the stress beyond which the sample begins flowing). 



3.2 Slow dynamics and aging in some jammed systems 

In their work on jammed fragile systems, Cates et al. consider hard spheres 
large enough that their thermal motion can be neglected [38]. Liu and Nagel 
[22], while proposing the jamming phase diagram for supercooled liquids, 
closely-packed grains and bubbles, speculate on the role of interparticle in- 
teractions in the jamming phenomenon. These predictions have been veri- 
fied experimentally by Trappe et al. who construct a jamming phase diagram 
(adapted in Fig. lb) for attractive colloidal particles [59]. The onset of jam- 
ming in these experiments is identified with the appearance of a low-frequency 
plateau of the elastic modulus that signifies the formation of an intercon- 
nected solid network [60] that is reminiscent of the stress-supporting force 
chains suggested in [38]. Owing to the attractive interparticle interactions 
that restrict the mobility of the individual constituents, jamming is found to 
set in at smaller values of the particle volume fraction than in the case 
of hard spheres. The jamming-unjamming transition in this system is there- 
fore strongly dependent on the interparticle interaction strength U. Like in the 
case of hard spheres, jammed attractive particles can also be unjammed by the 
application of stresses or by decreasing <fi. The composite phase diagram pro- 
posed in [59] has a suitably normalised applied stress a/a (er = k B T/a 3 , 
where a is the radius of the colloidal particle) and k B T/U as the three axes, 
such that an appropriate increase in one or more of these parameters results 
in fluidisation/ unjamming. 

A feature of systems near the jamming transition, such as supercooled liquids 
and glasses, is the presence of dynamical heterogeneities [20,61]. The slow 
and heterogeneous dynamics in jammed systems such as sand [62], colloidal 
hard sphere glasses [63] and coarsening foams [64] have been studied in de- 
tail. Due to the broad distribution in the depths of the energy barriers in the 
potential energy landscape characterising these systems, the mechanical prop- 
erties of SGMs evolve or 'age' continuously. Soft colloidal pastes formed from 
jammed microgel particles show many history-dependent properties that are 
interpreted in terms of aging phenomena [6,42]. Their strain recovery depends 
on the wait time (age) t a after flow cessation and can be scaled onto a master 
curve for all applied stresses below the yield stress a y . This type of scaling is 
characteristic of aging samples. This aging behaviour can be interrupted by 
applying large enough shears, a process called rejuvenation. Perhaps the most 
important result of this work is the observation that the amplitude of the stress 
uniquely determines the long-time memory and slow evolution of the colloidal 
paste. Microrheological measurements of the time-dependent creep compliance 
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of an aging cytoskeleton [26], which essentially comprises a crowded network 
of semi-flexible biopolymers, can be scaled onto a master curve very similar 
to the kind observed for microgel pastes [6]. The aging behaviour, intermit- 
tency and slow dynamics observed in this living system strongly resemble the 
dynamics observed in inert soft glasses. 
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Fig. 3. XPCS data (denoted by circles) for the intensity correlation function g2(q, t) 
for 3.0 wt.% laponite at t a = 90000s measured at q = 0.135 nm^ 1 is shown in a). The 
solid line is the result of a fit to the form 52(9,*) = 1 + /3[exp(— (i/T s ) L5 )] 2 , where 
(3 = bA 2 depends on the Siegert factor b and the amplitude of the short-time motion 
through A as discussed in detail in [24] . The dashed line in a) shows a poor fit of the 
data to a simple exponential and is meant to highlight the compressed exponential 
nature of the observed decay of 32 (9, t). The q dependence of the relaxation time 
t s at t a = 40000s is described by the relation t s ~ g _1 and is shown in the inset of 
a). Part b) shows a plot of log(r s ) vs. log(t a ) at q = 0.135nm _1 . Superlinear power 
law aging (r s ~ t^ 8 ) over the entire range of accesible t a s is observed. Some of the 
data for this figure has been adapted from [24]. 

Soft materials at the jamming transition often exhibit aging phenomena and 
universal, non-diffusive relaxation processes [2,20]. Interestingly, for systems 
as diverse as fractal colloidal gels, concentrated emulsions, micellar polycrys- 
tals, lamellar gels and synthetic clay suspensions [2,3,7,8,24], the DSF f(q,t) 
measured by PCS exhibits a two step decay. In spite of the arrested dynam- 
ics expected in these systems, f(q,t) is found to decay completely. In some 
cases, the characteristic relaxation time Ty measured for the faster, initial de- 
cay scales with q as 17 ~ q~ 2 , indicating caged diffusive dynamics [2]. In the 
case of a fractal colloidal gel, the DSF f(q,t) characterising the initial decay 
exhibits a stretched exponential form with an exponent of 0.7 [65] and is un- 
derstood in terms of the contribution of thermally excited, overdamped modes 
of the gel strands. The second, slower relaxation in these soft glassy materials 
is more counter-intuitive and is characterised by a DSF with a compressed 
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exponential form (f(q,t) ~ exp[— (t/r s ) p ], where p ps 1.5), with t s varying 
with q according to the relation t s ~ q^ 1 . As an example, Fig. 3a shows the 
intensity autocorrelation function g2(q,t) (indicated by circles) measured at q 
= 0.135 Jim -1 and t a = 90000 s for a 3 wt.% laponite clay suspension [24]. 
The solid line shows a fit to g2(q,t) corresponding to the following form for 
f(q,t): f(q,t) ~ exp[— (t/r s ) 15 ] [24]. This compressed exponential form of the 
DSF describes the data remarkably well when compared to a fit of the same 
data to a simple exponential form (indicated by the dashed line). The charac- 
teristic time scales r s when plotted vs. q follow the relation r s ~ q~ l at every 
sample age (inset of Fig. 3a shows the data at t a = 40000s). Evidence of such 
non-diffusive slow relaxation processes in aging laponite suspensions, but at 
the much smaller wave vectors accessible by DLS, have also been reported in 
[8]. 

In sharp contrast to the slow relaxation processes in supercooled liquids whose 
DSFs are characterised by stretched exponential forms f(q, t) ~ exp[— (t/r s ) p ], 
where p < 1 and r ~ 1/q 2 (with corrections for de Gennes narrowing [66]), the 
observed faster-than-exponential relaxation of these soft glassy materials is in- 
consistent with diffusive relaxation processes. Instead, as described in Section 
2.1, this unusual form of the DSF and the inverse linear dependence of r s on 
q are understood in terms the ballistic motion of elastic strain deformations 
that occurs in response to slowly developing localized sites of stresses [3,40]. 
Indeed, the picture developed heuristically by Cipelletti et al. [2] and devel- 
oped by Bouchaud and Pitard [40], in which these sites are identified with 
randomly positioned stress dipoles whose intensity grows linearly in time, ac- 
counts accurately for the form of the DSF and its wave vector dependence. 
However, the large array of materials in which these dynamics have observed 
experimentally [2,3,7,8,24,67] points to a more general principle underlying 
the formation of such sites of stress in jammed systems and calls for further 
investigation, both theoretically and experimentally, into this phenomenon. 

A common feature of these non-diffusive slow dynamics is the steady evolu- 
tion, or aging, of the characteristic relaxation time r s . However, in contrast 
to the compressed exponential shape of f(q,t), which is remarkably univer- 
sal, the precise nature of the aging behavior is strongly system dependent. 
For instance, t s varies as a power law with system age t a for colloidal frac- 
tal gels, micellar polycrystals, and lamellar gels, with power-law exponents in 
the range between 0.4 and 0.9 [2], in rough agreement with the prediction of 
the Bouchaud-Pitard model, r s ~ t 2 J z . However, experiments on compressed 
emulsions [2] show a superlinear dependence of r s on t a , in sharp contrast to 
the predicted behavior. As shown in Fig. 3b, the XPCS studies on laponite 
clay suspensions [24] similarly display superlinear behavior with r s ~ t^ 8 , 
where t a is the time since formation of the suspensions. DWS and XPCS mea- 
surements of aging depletion gels formed from silica nanoparticles also exhibit 
compressed exponential relaxation processes that evolve with time after the 
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Fig. 4. (a) Intensity autocorrelation function g2(q, t) at q = 0.14 nm _1 for a 3 wt. % 
laponite suspension initially aged for 3 days and then rejuvenated by filtering. Re- 
sults are shown for three times since rejuvenation: t a = 9600 s (circles), 28000 s 
(triangles) and 115000 s (squares). Solid lines are the results of fits to a compressed 
exponential lineshape with compression exponent 1.5. (b) Circles indicate the char- 
acteristic relaxation times extracted from g2(q,t) at q = 0.14 am" 1 as a function 
of time since rejuvenation. The dashed line has a slope of 1.0, indicating a linear 
relationship, t s ~ t a . 



cessation of strong shear as r s ~ t a [67]. In dynamic light scattering (DLS) 
studies of laponite and colloidal gels, an initial exponential increase of r s with 
t a is observed, followed by a regime of 'full aging' where t s varies as a power-law 
with t a [3,7,8]. Thus, no clear pattern has emerged for the temporal evolution 
of these non-diffusive dynamics. 



The sensitivity of the aging behavior of these slow dynamics to details of the 
system is further illustrated by recent XPCS measurements by our group on 
laponite suspensions following rejuvenation. Rejuvenation was accomplished 
by subjecting suspensions that have aged for long periods (3 days in this case) 
to strong shear flow by filtering them through 11 /im diameter pores. The sus- 
pensions initially re-fluidise due to the filtering, after which aging resumes and 
they progressively regain their solid-like consistency. Fig. 4a displays g2(q, t) at 
q = 0.14 nm _1 for a 3 wt.% laponite suspension measured at three times since 
rejuvenation. The solid lines are results of fits to the compressed exponential 
form with exponent 1.5, identical to the line shape observed for laponite sus- 
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pensions during aging after their initial formation, as shown in Fig. 3a. The 
characteristic relaxation times t s at q = 0.14 nm" 1 extracted at various times 
since rejuvenation t a are shown in Fig. 4b. The dashed line in the figure has 
slope 1.0, indicating t s ~ t a , clearly distinct from the XPCS results for initially 
formed laponite, t s ~ t^ 8 , shown in Fig. 3b [24]. We attribute this contrasting 
aging behaviour to differing origins of local stress in each case. As described 
in [24], sites of developing stress in the initially formed suspensions are iden- 
tified with a growing interparticle repulsion due to charge dissociation. In the 
rejuvenation, on the other hand, we identify them with residual local stress 
that is loaded into the sample by the strong shear flow and that is relieved 
by slow particle rearrangments. This more heterogeneous stress distribution, 
while leading to the same distribution of strain velocities and hence to the 
same form of f(q,t), evidently evolves differently in time, creating distinct 
aging behavior. 

We note, however, that these new XPCS results on laponite rejuvenation sug- 
gest a need to re-evaluate the role of growing interparticle repulsions in the 
formation of soft solid laponite suspensions. As described in [24], a key piece 
of evidence for such growth was a steady increase with age of the short-time 
plateau value of g2(q, t) and the increasingly constrained caged particle motion 
inferred from this increase. As Fig. 4a illustrates, laponite suspensions follow- 
ing rejuvenation unexpectedly display a very similar increase in the plateau 
value despite the fact that the interparticle interactions are likely unaffected by 
the rejuvenation. One possibility for this increase in the rejuvenated laponite 
suspensions is a slow strengthening of force networks in the suspensions weak- 
ened by the shear flow and a concomitant decrease in thermal fluctuations 
of these networks that create a fast partial decay in g2{q,t) in analogy with 
the picture developed in [65]. Thus, the increase in the plateau value during 
formation and after rejuvenation could have different origins. However, the 
strong similarity of the two behaviors suggests otherwise, calling into question 
the interpretation introduced in [24]. Further study of laponite rejuvenation 
with XPCS would be helpful in resolving this issue. Nevertheless, the contrast 
between Figs. 3b and 4b illustrates clearly how the aging behavior of the slow, 
non-diffusive dynamics is sensitive to the precise source of the local stress. 



3.3 The SGR model: experimental status 

Glassy rheology is a common feature in soft materials characterised by metasta- 
bility and disorder. A signature feature of soft glassy rheology is the weak 
power law behaviours of the viscoelastic moduli. Experimental systems that 
verify the power law behaviour of soft glassy materials include concentrated 
microgel solutions [31], foams [32], paint [33] and compressed emulsions [34], 
all of which have G's that show a weak power law frequency dependence, with 
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the exponent (x-1) of the power law lying in the range 0.1-0.3. That is, the 
noise temperature x within the SGR model lies in the range 1.1-1.3. In the 
case of the compressed emulsions [34], G' is always significantly greater than 
G" above the random close packing fraction (4> rcp ~ 0.635) of the droplets. As 
the droplet volume fraction is increased beyond rcp , the droplets squeeze 
against each other and deform, and the resulting constraints on the droplet 
dynamics manifest as an increasingly flat elastic response. The elastic moduli 
of hard sphere colloidal suspensions, over a limited range of frequencies, also 
show power law frequency dependences characterised by exponents that de- 
crease monotonically as the colloidal volume fraction increases toward the 
kinetic glass transition <p g = 0.58 [68]. These results are understood in terms 
of the kinetic arrest of particles trapped in their cages as the kinetic glass 
transition is approached. This decrease in the power law exponent towards 
(x — 1) =0 indicates an approach towards the glass transition of the SGR 
model and thus directly connects that transition, at least qualitatively, to the 
kinetic glass transition in hard sphere colloids. Further, for volume fractions 
above <f> g , colloidal suspensions display aging [69]. 

Creep experiments on onion phases (20% AOT in a 15 g/1 brine solution) also 
uncover very slow relaxation processes [70]. The frequency dependent moduli, 
which track one another over more than three decades of angular frequency, 
show two distinct regimes: a low-frequency power law regime with an exponent 
close to zero that describes the slow relaxation of grain boundaries and long 
length-scale reorganisation processes, and a high-frequency regime which is 
dominated by the response of individual onions. The first part of the response 
is strongly reminiscent of SGR. The elastic and viscous moduli of suspensions 
of aging laponite, after showing Maxwellian behaviour at the earliest times, 
eventually exhibit power law frequency dependences, with the exponent of the 
power law becoming smaller as the sample ages [71]. This observation of aging 
in laponite when exponent (x — 1) > at first glance seems at odds with 
the SGR picture, since x > 1 constitutes the ergodic regime in the model. 
However, as mentioned above, studies of laponite indicate that the aging is 
driven at least partially by an evolving interparticle potential [24], a feature 
that clearly lies outside the range of the model. 

Another system exhibiting signature features of soft glassy rheology is the 
thermotropic smectic liquid crystal octylcyanobiphenyl (8CB) confined within 
the pores of a highly porous and compliant hydrophilic aerosil gel (size of in- 
dividual aerosil particles is 7 nm) [35]. In this system G' and G" show power 
law frequency dependencies that arise from an interplay of the smectic order 
and the quenched disorder introduced by the gel. For example, Fig. 5a dis- 
plays G'(uj) at 21°C (red circles), several degrees below the nematic to smectic 
transition, and at 36°C (blue squares), above the transition. At high temper- 
ature, G' is dominated by the colloidal gel. The greatly enhanced modulus in 
the smectic phase at lower temperature includes a weak power law compo- 
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nent whose exponent (x — 1) decreases with increasing aerosil density p s , as 
shown in the inset to Fig. 5a. This decrease thus signals an approach to the 
glass transition of the SGR model, x — 1, [29,30] with increasing quenched 
disorder (quantified by p s ). The key role of the smectic order in this rheology 
is illustrated by contrasting this behavior with the rheology of the simple or- 
ganic liquid dibutyl phthalate confined to the pores of an aerosil gel, shown 
in Fig 5b. G'(u) measured at 21°C (red circles) and 36°C (blue crosses) are 
essentially identical and show no frequency dependence as expected for the 
rheology of a colloidal gel in a low viscosity solvent. Scaling behavior of the 
smectic modulus for the 8CB+aerosil indicates that the enhanced elasticity 
and glassy rheology result from a dense network of screw dislocations that 
arise due to the quenched disorder and whose motion is pinned by the dis- 
order [35]. An excellent overlap of the data for different disorder strengths is 
observed when the ratio ^ X ^ G „ is plotted vs. temperature below the pseud- 
ocritical region [35] . These observations point to a one-to-one correspondence 
between p s and x, where x quantifies the noise temperature influencing the 
dynamics. 

Microrheological measurements in living cells show power law frequency de- 
pendences of the viscoelastic moduli of the cytoskeleton, where the effective 
noise temperature lies in the range 1.1 < x < 1.3 and is found to depend 
strongly on the nature of the biological intervention employed in the experi- 
ment (such as the addition of a contractile histamine) [27]. Those biological 
processes that cause contractile acivity in the cell or result in the polymerisa- 
tion of the cell proteins can drive the system towards a glass transition. Dis- 
order and metastability, therefore, contribute significantly to the mechanical 
functions of the living cell, with the noise temperature x serving as a useful 
measure of its deformation and flow behaviour. Similar glassy behaviour is 
also observed in the rat airway smooth muscle during actin modulation [72]. 
The scale-free soft glassy rheology of living cells has been a subject of recent 
research [73,74,75]. 

Under continuous shear flow, the potential energy landscape of soft glassy 
materials gets modified, such that the depths of the minima decrease tem- 
porally according to the relation E(t) = E — ^k^t 2 , where k is the elastic 
constant relating the stress to the deformation and 7 is the macroscopic strain 
rate. Many SGMs, such as concentrated emulsions [76] and foams [77] show 
Herschel-Bulkley type flow curves (discussed in section 2.2) with a well-defined 
yield stress a y . For stresses of magnitude a > cr y , the flow curve assumes a 
power law form with a power (x-1) of around 0.1. Another example is a micro- 
gel paste made of jammed polyelectrolyte particles that can flow when sheared 
above a certain yield stress value. The glassiness of this material is evident 
from its Herschel-Bulkley flow curve [42]. The synthetic clay laponite [71] also 
shows a power law flow curve with noise temperature x = 1.20, a value slightly 
inconsistent with the values of x obtained from the frequency-dependent mea- 
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Fig. 5. Rheology data for the elastic moduli G' of 8CB dispersed in a porous aerosil 
gel of density p s = 0.10 g/cm 3 is shown in a). The red circles correspond to G' 
measured in the smectic phase (T = 21°C) which shows a power law contribution 
characterised by an exponent (x — 1). The solid line through the data at 21°C is the 
result of a fit to a weak power law component plus a frequency-independent com- 
ponent as detailed in [35] . The viscous modulus G" in the smectic phase also shows 
an identical power law dependence when the contribution of the gel is subtracted 
from the response (not shown in this figure). The blue squares correspond to the 
approximately frequency independent behaviour when 8CB is in the nematic phase 
(T = 36°C). The plot of (x — 1) vs. p s is shown in the inset. The errors bars are 
computed from the temperature variation of {x — 1) at each p s . This data is dapted 
from [35]. Part b) shows the frequency response of a simple organic liquid dibutyl 
phthalate confined to the pores of an aerosil gel (p s ~ 0.09 g/cm 3 ). The G' mea- 
sured at 21°C (red circles) and 36°C (blue crosses) are both frequency independent, 
while their magnitudes are almost identical to the G' measured for the 8CB-aerosil 
mixture (with p s ~ 0.10 gm/cm 3 ) when 8CB is in the nematic or isotropic phase. 



surements of G' and G" . This discrepancy leads the authors to conclude that 
the nonlinear flow behaviour of colloidal laponite is significantly more complex 
than the predictions of the SGR theory [29,30]. There is, therefore, a need for 
more experimental studies to verify the predictions of the SGR theory for the 
nonlinear response of SGMs. 
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4 Concluding remarks 



In this paper, we have reviewed theoretical models used extensively to un- 
derstand the slow dynamics of soft and living matter, focussing in particular 
on two key concepts: jamming and soft glassy rheology. The main purpose 
of the paper has been to provide a overview of these concepts and to cata- 
logue some of the extensive experimental studies on a diverse range of soft 
and glassy materials that relate to these models. These concepts continue to 
undergo theoretical development, as exemplified by the recent work relating 
the zero-temperature jamming transition to k-core percolation [78] and the 
introduction of tensorial constitutive models that generalise SGR [79]. On the 
experimental side, the emergence of techniques such as multispeckle correla- 
tion spectroscopy and the availability of high-precision commercial rheome- 
ters have greatly facilitated the study of soft glassy materials. Simultaneous 
rheometry and PCS studies that can correlate the bulk flow behaviour of 
these materials to their microscopic structure and dynamics are a potentially 
fruitful avenue for future experimental work. Such continued innovation in ex- 
perimental methods should lead to a better understanding of the aging, slow 
relaxational behavior, and glassy rheology observed in soft and living systems. 
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